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The purpose of this study is to describe a perfect fluid matter distribution that leads to a constant 
curvature region, thanks to the effect of a non-minimal coupling. This distribution exhibits a 
density profile within the range found in the interstellar medium and an adequate matching of 
the metric components at its boundary. By identifying this constant curvature with the value 
of the cosmological constant, and superimposing the spherical distributions arising from different 
matter sources throughout the universe, one is able to mimic a large-scale homogeneous cosmological 
constant solution. 
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I. INTRODUCTION 

Contemporary cosmology faces currently major chal- 
lenges. These include, for instance, the need to firmly es- 
tablish the existence and nature of dark energy and dark 
matter, or to find evidence of a more encompassing the- 
ory of gravity. The nature of dark energy [1] has been dis- 
cussed in several proposals, including quintessence mod- 
els, -which consider the slo-w-roll of a scalar field [2, 3], av- 
eraging of inhomogeneities at a cosmological scale yield- 
ing an effective scalar field [4], amongst other candidates. 
A possible unification of both "dark" components has 
also been suggested, resorting to a scalar field model [5] 
or an exotic equation of state, as featured by the so-called 
generalized Chaplygin gas [6] . 

Another outstanding difficulty in the context of the 
General Relativity (GR) description of the universe dy- 
namics and its relationship -with the standard model of 
particle physics is the cosmological constant problem, a 
notorious adjustment issue that so far lacks a proper ex- 
planation (see e.g. [7] and references therein for a dis- 
cussion and a list of the most -well kno-wn proposals for 
a solution). A different approach assumes that no extra 
energy content is needed, but that the fundamental la-ws 
and tenets of gravitation are incomplete, such that -with 
current kno-wledge just offers a lo-w-energy perspective of 
a more fundamental theory; as a consequence, modifi- 
cations of the Friedmann equation to include higher or- 
der terms in the energy density p (see e.g. [8, 9] and 
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references therein) have been proposed or, at a more 
fundamental level, changes to the action functional. A 
rather straightforward approach lies in replacing the lin- 
ear seal ar curvature term in the Einstein-Hilbert ac- 
tion -with a function of the scalar curvature, f{R). Al- 
ternatively, one could resort to other scalar invariants 
of the theory [10]: extensions relying on a functional 
dependence of the action on the Gauss-Bonnet invari- 
ant G ^ H^ - iRf.^Rt^" + Rapf.uR"^'"' are perhaps the 
most -well studied models, as this could arise in a lo-w- 
energy effective description of String Theory and carry 
several implications in brane-world scenarios [11]. The 
more tractable class of f{R) theories have had consider- 
able success in replicating the early period of rapid ex- 
pansion of the universe, -with the Starobinsky infiation- 
ary model f{R) = R-\- aR^ being the most -well- kno-wn 
proposal [12] . The accelerated expansion of the universe 
has also been the subject of attention [13]. Solar system 
tests could also bring further insight, mostly arising from 
the parameterized post-Ne-wtonian (PPN) metric coeffi- 
cien ts derived from this extension of GR. Ho-wever, some 
disagreement exists in the community, with some argu- 
ing that no changes are predicted at a post-Newtonian 
level (see e.g. [14] and references therein); amongst other 
considerations, this mostly stems from an approach based 
either in the more usual metric affine connection (that is, 
where the affine connection is taken a priori as depend- 
ing on the metric), or in the so-called Palatini approach 
[15] (where both the metric and the affine connection are 
taken as independent variables). As an example of a clear 
phenomenological consequence of this extension of GR, it 
has been shown that f{R) = /o-R" theories yield a gravi- 
tational potential which displays an increasing, repulsive 
contribution, added to the Newtonian term [16]. 

Notwithstanding the significant literature on these 
f{R) models, steps have been taken to explore another 



interesting possibility: not only that the curvature is 
non-trivial in the Einstein-Hilbert Lagrangian, but also 
that the coupling between matter and geometry is non- 
minimal, i.e. enforced solely by the invariant ^/—gd'^x, 
the use of the metric to raise and lower indexes and the 
associated covariant derivative; indeed, in GR covariantly 
invariant terms in Cm should be constructed by contrac- 
tion with the metric {e.g. the kinetic term of a real scalar 
field, g^^4'^^4'^v)'- a non-minimum coupling would imply 
that geometric quantities (such as the scalar invariants) 
would explicitly appear in the action. This possibility 
can have deep phenomenological implications: in par- 
ticular, it would imply that regions with high curvature 
(which, in GR, are related to regions of high energy den- 
sity or pressure) could lead to considera ble deviations 
from the dynamics predicted by Einstein's theory [17]. A 
wide range of results has arisen from this hypothesis, in- 
cluding the formal equivalence with a multi-scalar-tensor 
theory [18], the impact on solar observables [19], the pos- 
sibility to account for dark matter [20] , post-inflationary 
reheating [21] or the current accelerated expansion of the 
universe [22]. 

When studying the latter, it was noted that one can- 
not recover an exponential expansion due to a constant 
curvature, i.e. a universe dominated by a cosmological 
constant (CC), when working under the assumption of 
a Friedmann-Robertson- Walker metric. Since this is the 
simplest of all culprits behind the mentioned accelerated 
expansion, one now reassesses this issue and attempts 
to obtain a local, spherically symmetric solution for the 
interstellar matter (ISM) distribution that, in the con- 
text of a non-minimally coupled model, gives rise to a 
constant scalar curvature. This yields a mechanism that 
generates what can be viewed as a local CC. The super- 
position of several spherical distributions throughout the 
universe could then lead to a universal CC. Of course, 
this is not a solution for the cosmological constant prob- 
lem, but a proposal to account for the observations once 
a mechanism is found to make all other contributions 
vanish. 

The goal of this work is then to achieve this descrip- 
tion without the requirement of a specific model for the 
non-minimal coupling or a somewhat awkward and con- 
voluted density profile — i.e. although one shows that 
only a quite specific distribution gives rise to a constant 
curvature, this can be regarded as an admissible profile 
for the ISM medium. One may also regard this work as 
yet another example of the strong impact a non-minimal 
coupling can have in the description of matter distribu- 
tion, when compared with the standard results arising 
from GR. 

Related proposals aiming to describe, for instance, 
dark energy through some variation of the concept of 
a non-minimal coupling between matter and geometry 
have been put forward previously, addressing the prob- 
lem of the accelerated expansion of the universe [23] and 
the existence of a CC [24]. However, one posits that the 
approach proposed here is somewhat more economical 



and elegant, as it does not explicitly introduce any addi- 
tional matter field, and it does not set any specific form 
for the non-minimally coupled terms. 



II. THE MODEL 

One considers a model that exhibits a non-minimal 
coupling between geometry and matter, as expressed in 
the action functional [17], 



S= f [kR + f2{R)C] ^^d* 



(1) 



Variation with respect to the action yields the modified 
Einstein field equations. 
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with K = c^ /IQttG, Af^i, = V^V,y — g^i^O and C = — p 
(see Ref. [25] for a discussion about this choice within the 
context of a non- minimal coupled model). As expected, 
GR is recovered by setting fi{R) = 2kR and f2{R) — 1- 
Resorting to the Bianchi identities, one concludes that 
the energy-momentum tensor of matter may not be con- 
served in a covariant way, as 
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Again, in the absence of a non-minimal coupling, 
/2 {R) = 1 and the covariant conservation of the energy- 
momentum tensor is recovered. This feature implies that 
the motion of the matter distribution described by a 
Lagrangian density C does not follow a geodesic curve. 
Clearly, a violation of the Equivalence Principle may 
emerge if the r.h.s. of the last equation varies signifi- 
cantly for different matter distributions, which suggests 
a method of testing the model and imposing constraints 
on the associated couplings. 

A constant curvature leads to i?^^ = i?o5/ii'/4 -^ R — 
Rq and L — — p; inserting this into Eq. (3), one gets 



{n + F2p)g^,Ro^4.F2A^,p-2f2T^, . (4) 

One writes the metric as 



ds^ =. -e^'^^-Uf + e2-('')dr2 + dfl^ 



(5) 



Assuming a static configuration, u^ = (uo,0); since 
u^,u'' = g"°ul = -1, one gets uq = (-500)^^^ = e"^'''' . 
Thus, the perfect fiuid form for the energy-momentum 
tensor 



while the (0, 9) component of the field Eq. (9) reads 



Tfj.„ = {p + p)u^,Uv + pg^i, -^T = 3p- p 
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reads 
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One assumes a constant equation of state (EOS) pa- 
rameter Q '^ uj — pI p < 1, where w = 0, 1/3 and 1 
gives dust, ultra-relativistic and ultra-stiff matter, re- 
spectively. Taking the trace, one obtains 



(K + F2p)i?o--3F2np+ 2/2(1 -3w)p . (8) 



Replacing into Eq. (4), one gets 
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Notice that both equations show that the only constant 
solution is p = — 7- -Ro = 0, as in GR (unless f2{Ro) = 
0, which is unphysical). If F2 ^ 0, equating the two 
equations above yields 
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+ a'] p' ^ p' = KRoe^r 
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where X is a constant with dimensions of density. Notice 
that K = leads top = Oorw = — 1, the GR conditions 
for constant curvature. 

Adding Eqs. (13) and (14) yields 
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One then has 
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/2(l+w)pe2- = 2F; 
which, again using Eq. (13), reads 
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implicitly defining p in terms of v. 

In order to obtain a single differential equation for p, 
one first isolates the metric coefficients as 



e^- = 



KRor 
p' r 

1 1+UJ 

r 



(18) 



/2 



1 



2 F2iKRor) 
l + u; /2 
2{KR^rY F2 



7PP 



ip'r+pp" 



-pp 

r 



One now takes the definition of the scalar curvature. 



so that 
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where p' = p^^ and one assumes only a radial dependency 
for all quantities. 

Since e^'^500 + e^^grr — 0, one obtains from Eq. (9) 
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and, replacing Eq. (18), obtains a differential equation 
for p. 
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In the case of the minimal coupling f2{R) = 1, this is 
ill-defined unless p' — 0, as discussed before. This is also 
attained by setting K = 0. 



A. Boundary matching conditions 

The purpose of this work is to show that one may lo- 
cally mimic a CC via a suitable matter distribution, in 
the context of a non-minimal coupling between the latter 
and curvature. In line with the Cosmological Principle, 
it was argued in the introduction that, even in a cosmo- 
logical context, A arises not as a fundamental constant of 
nature, but from an averaging of all the surrounding mat- 
ter distributions with constant curvature (perhaps even 
each with a different curvature) . A more careful study of 
the averaging of matter distributions with different val- 
ues of the (constant) curvature is delayed for a future 
work (see e.g. [26], as well as [27] for a discussion of 
the Einstein-Strauss and Lindquist- Wheeler solutions for 
the issue of embedding a local "bubble" in an expanding 
spacetime issue in GR). 

However, one cannot simply assume that such averag- 
ing yields a constant density space surrounding the in- 
terior solution under scrutiny: from Eqs. (13) and (14), 
one sees that the only constant solution is p — Q: this is 
expected, as this work shows that a constant curvature 
enforces a specific non-flat profile for the density. 

In order to circumvent this issue, one assumes that the 
interior solution is surrounded instead by an exterior vac- 
uum solution with a constant curvature R = 4A. In this 
region, p = and the effect of the non-minimal coupling 
vanishes: one recovers the usual Schwarzschild-de Sitter 
metric. 
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so that 



Inserting p"{rs) = into Eq. (23), one obtains a relation 
between the total mass M and r^ , independent of A and 
other parameters: 
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Notice that M is not necessarily the gravitational mass, 
but merely an integration constant of the field equations 
in vacuum. 

Using Eq. (15) to evaluate Eq. (20) at r = r^, one 
obtains 
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Continuity of v is imposed by an appropriate integra- 
tion constant, since only its derivatives show in the dy- 
namics. Regarding a, matching at r = r^ shows that the 
constant K is related with the slope of the density at the 
boundary of the spherical body. 
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which is also obtained by evaluating Eq. (15) at r — r^ 
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In order to match the interior and exterior solutions 
at the radius of the spherical body r ^ Vs, defined as 
p{rs) =0, Eq. (18) lead to 
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III. NUMERICAL ANALYSIS 

In order to pursue a numerical analysis of the suggested 
scenario, one introduces the dimensionless scaling 
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thus implying that p"{rs) = 0; this does not depend upon 
the specific form of v and ct, but only on the boundary 
conditions i^'irg) ~ —a'(rs) and p{rs) = 0. Indeed, in- 
serting these into Eq. (13) leads to 



having used 
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with dn — c/Hq — 4.22 Gpc the Hubble distance. 

Notice that a negative sign is used in the first definition 
of Eq. (30), so that a decreasing density p' <Q ^f K <Q 
[cf. Eq. (15)) yields a positive 6{x) function. With this 
transformation, Eq. (20) only involves the parameters e 
and A, 
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where the primes now indicate differentiation with re- 
spect to X. 
The boundary conditions for the density translate into 



0(1) =0 
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(34) 



as can be checked by evaluating Eq. (33) at a; = 1 and 
taking e"{l) = p"{rs) = 0. 

The numerical procedure is simple: for a given K, Rq, 
u) and /2, compute e and A. Evolve the above differential 
equation between x ~ and a; = 1, with the boundary 
conditions conditions above, and then read the density 
from Eq. (30). Notice that the constant K does not ap- 
pear in Eq. (33); one may view it as a scaling parameter 
set by the central density of the spherical matter distri- 
bution, 



Pe = --K9{0) 
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l + 3A-e 
since one expects the density to drop and K < 0. 

A. Smooth transition 

Before continuing, one addresses another possibil- 
ity: instead of imposing continuity with an outer 
Schwarzschild-de Sitter metric, one could perhaps define 
the boundary of the spherical mass distribution as a ra- 
dius r = rout where p{rout) = Pout and p'{rout) = 0, i.e. 
where a smooth transition to an exterior solution occurs. 
Notice that rout 7^ Ts, the latter being defined through 
Eq. (28) and characterizing the interplay between the 
model parameters. However, this cannot be done, as it 
can be seen from Eq. (15) that p'{rout) = — >■ rout) = 
or K ^ Q] the latter yields a constant p — pout sphere 
with no physically meaningful boundary with the iden- 
tical distribution surrounding it, and which only allows 
for a vanishing constant curvature i? = 0, as discussed 
before. 



B. Metric elements 



One may rewrite Eq. (18) in terms of the redefined 
quantities given by Eq. (30), 
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As discussed, these do not depend upon the integration 
constant K, showing that the direct physical meaning of 
6 and its derivative is in the metric itself. 

As discussed above, the boundary conditions and their 
scaling may be extracted from the matching with the 
rescaled Schwarzschild-de Sitter metric, after substitut- 
ing M = (167r/3)Krs, which reads 
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By matching the two sets of equations, one again recovers 
the boundary conditions Eq. (34), as expected. 
From Eq. (37), one has 
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Integrating the first of Eq. (36), one finds that 
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setting the integration constant that provides the match- 
ing of this metric clement at the boundary r = rg , as was 
discussed following Eq. (28). Assuming a smooth 0{x), 
one can see that no singularity arises, since inspection 
of the above expressions shows that neither e^°' nor 
vanish. 
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C. Outer Schwarzschild-de Sitter metric 

As discussed before, the mass parameter in the outer 
Schwarzschild-de Sitter metric does not need to refiect 
the gravitational mass of the perfect fluid matter distri- 
bution, as it is an integration constant of the Einstein 
equations in vacuum. As a result, one has shown that 
continuity of the inner and outer metric requires that 
M = (167r/3)Kr,. 



The gravitational mass Mg is defined as 



Mg = 47r / pr dr = Airpc 
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If one aims to circumvent this distinction between M 
and the total mass of the spherical body, the condition 
Mg = M yields 



with 
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Resorting to Eq. (31), one finds that this sets the 
central density of the spherical body as 
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IV. MIMICKING A COSMOLOGICAL 
CONSTANT 

One now searches for a solution that mimics a cosmo- 
logical constant from a typical dust distribution (w « 0), 
i.e. one sets i?o = 4A. From Eq. (30), this translates 
into the simple relationship A = e. 

The boundary of the spherical body corresponds to 
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A<-l/2 . 

The condition r^ <C dn becomes 
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-1 0.57 Gpc 

-2 1.97 Gpc 

-10 19.8 kpc 

-20 0.20 pc 

-30 0.41 AU 

-40 612 km 

-50 6.12 m 



TABLE I: Table with radius of spherical distribution for dif- 
ferent values of A = e ^ 0. 
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All obtained solutions exhibit a smooth transition from 
|0a(O)| ~ 1 to 0a(1) = 0, and 0^(0) = 0. This good 
behaviour is illustrated by the central value, given by 
6'a(0) « 0.9 for all A = e < 0.1. This similarity enables 
one to fit them through 
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Notice, however, that the above fitting is merely illustra- 
tive, as it does not fulfill the required boundary condi- 
tions. 

One may now set the constant K , 
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Since Eq. (44) may be approximated by 
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= and 9' {I) = - 


\/6. Its 
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with boundary conditions 9{1) 
(normalized) solution is plotted in Fig. 1. 

Numerically, one finds that the solution for Eq. (33) 
with different values of A = e ~ (dubbed 9x{x)) are 
all very similar to the solution ^0(2;) of Eq. (46) (where 
A = e = 0); this is visible in the inset of Fig. 1, which 
plots the normalized ratio between 9\ and 9q, 
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one may rewrite the above as 



X«-0.19(^) p. 
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The first condition of Eq. (51) is illustrated in Table I. 



A. M = Mg condition 

Recalling that the condition M — Mg implies a spe- 
cific radius for the dust distribution, one writes 




FIG. 1: Main: scaled solution 6lo(x)/6lo(0) for e = A = 0. 
Inset: ratio Qx for e = 10"'^ (full), 10"^ (dashed) and 10"^ 
(dotted). 




FIG. 2: Comparison of the metric coefficients e^" (full) and 
e^'^ (dashed) derived from the e = solution (black) and the 
usual GR interior solution (grey). 






-11.4 Mpc 



(53) 



having used the solution 9{x) with e = A = (i.e. 
^0(2;) instead of ^^(2;)), for simplicity, in the numerical 
integration of /, and pn = 1 H atom/cm'^ = 1.660 x 
10^21 kg/m^. 

The typical distance between galaxies is of the order 
of Ts ~ IMpc {e.g. Andromeda is 778 Kpc away); one 
thus obtains 



130 H atom/cm^ 



(54) 



close to the lower bound of HII star-forming re- 
gions. The interstellar medium attains densities up to 
10** H atom/cm^, allowing for a radius of the constant 
curvature bubble down to ^ 100 Kpc. 

In the A = e ~ regime, Eq. (43) becomes 
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yielding a lower bound for pc well below the typical min- 
imum density of the interstellar medium. 



B. Metric elements 



With A = e ^ 0, the metric elements read 



One may compare this with the usual metric derived 
from GR for an interior solution, assuming for simplicity 
a vanishing pressure dust distribution, p = 0: 



771(3 (r) — 4:TT j p{r)r^dr , 


^2. ^ ,^_ma{r^^-' 



(57) 



STTKr 



exp 



mair) 



Q r (STTKr — mcir)) 



dr 



so that M ~ Mq — mc{rs). In terms of the rescaled 
quantities given by Eq. (30), one reads 
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having used M = (167r/3)Krs. 

It is advantageous to evaluate analytically the above 
integrals; for this, one resorts to the quartic cubic fit 
of the e = A — solution Eq. (48). One finds that 
the obtained metric coefficients are rather different from 
their GR counterparts (matching only at a; = 1, due to 
the condition M = mg{x = 1). However, this is to be 
expected: our solution arises from R = Rq = const., 
while in GR one has R = —T = p. 



o2<T(a;) 



^2^{x) 
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I9'{x)] 
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x 




x2 

yexp 


2 



—e{x)9'{x)dx 



(56) V. SPECIFIC MODELS 

One now looks at specific models that enable the mim- 
icking of a cosmological constant, i.e. 



Z. Quadratic coupling 



A = 



F2A 



1+UJ /2 



(59) 



Instead of obtaining the value for e, one discusses the 
range of model parameters derived from the condition 
e ^ for several forms of J2{R)- One sets w = 0, typical 
of a pressureless dust distribution; however, the above 
shows that the effect of a non-vanishing EOS parameter 
does not change the order of magnitude of the results. 



A. Power-law^ coupling 



One first assumes 



/.m^i + l|- 



so that 



(60) 



One now aims at obtaining spherical mass distributions 
of constant curvature arising from a quadratic correction 
n = 2; one imposes that this correction is below the linear 
one previously addressed at the inflationary stage, when 
R^ Ri < [Mpc/hf, 



^\ «- ^3L X ini2^ 

i?2 > ^=Ml X 10"*^ ^ X 10^2 m.-^ 



From Eq. (59), one has 



(64) 



'-"'i 



576 



n^ 



dffR2 



< 2.3^x10 



-227 



, (65) 



showing that any subdominant quadratic coupling at in- 
flation allows for the local mimicking of the cosmological 
constant. 



F2R0 

V 



2n 



Rn 

Ro 



/to. 



(61) 



If n > 0, one obtains Rn S> Ro] negative exponents 
imply i?„ <C i?o- In the following, one explores specific 
values for the exponent n. 



3. Inverse power-law coupling 

Taking n = — 1, dark matter constraints require that 
i?_i - (10 Gpcy^ [20]. Thus, one obtains from Eq. (59) 



-Rn 



e^lT^y""^'""-'^"' 



Analogously, setting n — —1/3 implies R-1/3 
(10^ Gpc)-^ [20], leading to 



1. Linear coupling 

Replacing n = 1, one must take into account the con- 
straint arising from the Starobinsky scenario for inflation 
and the ensuing reheating process due to the presence of 
a linear coupling between curvature and matter, as de- 
scribed in Ref. [21], 






(62) 



with 1 < ^ < lO"'. Since i?i > i?o = 4A, and recaUing 
Eq. (32), one gets 



2 f R-1/3 

3 ^ Ro 
2 1 



1/3 



3(12r!A)'/' VlO^Gpc 



IH 



2 ( R-i/3dl 

3 V 12r2A 

2/3 



(67) 



8.6 X 10 



-5 



Both values for e are vanishingly small, showing that 
a inverse power-law coupling between matter and curva- 
ture can account for the rotation curves of galaxies and 
enable the description of spherical bodies with constant 
curvature Ro = 4A. 



B. Exponential coupling 



A 24r2A ins . ^ 



i?i dj^Ri 



One now assumes a power-law exponential coupling of 
the form 



which is vanishingly small, as intended, for any value of 
^ in the considered range. 



f2{R) =exp 



R 

Rn 



(68) 



so that 






2n 



Rn 



(69) 



This expression is identical to Eq. (61), and the conclu- 
sions apply, namely that either n ^ 0, i?„ >• i?o (positive 
n) or i?„ ^ Rq (negative n). 



VI. CONCLUSIONS AND OUTLOOK 

In this work one shows that the rich phenomenology of 
non-minimally curvature-matter coupled models allows 



for a striking result, namely the possibility of inhomoge- 
nenous spherical mass distributions with constant curva- 
ture. The density profile of such matter distribution was 
obtained under general assumptions, and shown to lead 
to physically admissible values for density and size, typ- 
ical of the interstellar environment. Finally, one estab- 
lishes the compatibility between the matter distribution 
here discussed and the previously obtained constraints 
on the parameters of a class of power-law or exponen- 
tial couplings /2(-R). Although the mimicking of a CC 
is, by construction, local, one tentatively argues that the 
averaging of several of the obtained matter distributions 
could lead to a universal CC, which hints a novel mech- 
anism for the mimicking of a CC type term. 
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